Electrostatic interaction of neutral semi-permeable membranes 
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We consider an osmotic equilibrium between bulk solutions of polyelectrolyte bounded by semiper- 
meable membranes and separated by a thin film of salt-free liquid. Although the membranes are 
neutral, the counter-ions of the polyelectrolyte molecules permeate into the gap and lead to a steric 
charge separation. This gives rise to a distance-dependent membrane potential, which translates into 
a repulsive electrostatic disjoining pressure. From the solution of the non-linear Poisson-Boltzmann 
equation we obtain the distribution of the potential and of ions. We then derive an explicit formula 
for the pressure exerted on the membranes and show that it deviates from the classical van't Hoff 
expression for the osmotic pressure. This difference is interpreted in terms of a repulsive electro- 
static disjoining pressure originating from the overlap of counterion clouds inside the gap. We also 
develop a simplified theory based on a linearized Poisson-Boltzmann approach. A comparison with 
simulation of a primitive model for the electrolyte is provided and does confirm the validity of the 
theoretical predictions Beyond the fundamental result that the neutral surfaces can repel, this mech- 
anism not only helps to control the adhesion and long-range interactions of living cells, bacteria, 
and vesicles, but also allows us to argue that electrostatic interactions should play enormous role in 
determining behavior and functions of systems bounded by semipermeable membranes. 
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I. INTRODUCTION 



It is hard to overestimate the role semipermeable mem- 
branes and osmotic equilibria associated with them play 
in our everyday life. The best known examples are the 
natural biological membranes, which are highly imperme- 
able to ions due to their phospholipid bilayer structure, 
but become semi-permeable when ion channels are open. 
Such lipid membranes with channel proteins surround 
all biological (eukaryotic and prokaryotic) cells [l[ . Many 
synthetic membranes used in electrochemical fuel cells 0| 
and dialysis Q take advantage of the semipermeability 
of their materials. The same concerns various types of 
synthetic vesicles , the promising gene and drug car- 
riers, as well as other systems used to mimic biological ob- 
jects. Being in contact with ionic solutions, such as gels, 
polyelectrolytes (including DNA, proteins, dendrimers), 
micelles, or colloids, a semipermeable membrane main- 
tains an unequal distribution of ionic solute concentra- 
tions, which leads to an ion density gradient across the 
membrane, determines its actual potential, and generates 
an osmotic pressure difference that the membrane has to 
sustain. This situation is traditionally referred to as a 
Donnan equilibrium 0, Q . 

A great deal of research has been devoted to under- 
standing the Donnan equilibria and the pressure exerted 
on semipermeable membranes. Although the problem 
was treated at various levels of sophistication, tradition- 
ally the reservoirs separated by a semipermeable wall are 
assumed to be sufficiently large, so that both solutions 
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FIG. 1: An electrolyte solution in equilibrium with a thin 
liquid film of thickness h bounded by neutral semiperme- 
able membranes. They are permeable for small counterions 
only (here anions), while the larger solute (here cations) are 
bounded to the region |a:| > h/2. The membranes have to 
support an excess osmotic pressure Ap(/i) which is the sum 
of the bulk Donnan pressure and an electrostatic disjoining 
pressure II(/i) in the interlayer. 



contain a phase with the bulk properties. Our paper 
deals with another situation, when two ionic electrolyte 
solutions consisting of large ions and small ions (re- 
ferred below to as counter-ions) are in equilibrium with 
a thin film bounded by the semipermeable membranes. 
As some mobile counter-ions will inevitably escape from 
electrolyte solutions, at some separation their clouds will 
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begin to overlap and give rise the effect similar to the cel- 
ebrated Derjaguin disjoining pressure [91 and will change 
dramatically the Donnan equilibrium [7| in the system. 
Some experimental observations [lo| support this idea. 
However, to our surprise, such a scenario, which repre- 
sents enormous interests for many biological and materi- 
als science problems associated with the membrane ad- 
hesion and long-range interactions and are also 
relevant to modern micro- and nanofluidics [13, [13 , has 
never been addressed before p^ . 

In our previous papers we made an attempt to calcu- 
late an excess pressure on a sernipermeable shell in con- 
tact with an inner 17| or outer [l8| solution of polyions. 
Since a non-linear Poisson-Bolzmann (NLPB) equation 
cannot be solved analytically for spherical geometry [l^ , 
its linearized version (LPB) has been used. Our present 
paper solves a pressure problem for a flat geometry of two 
interacting semipermeable membranes. We first solve 
semi-analytically a NLPB equation to evaluate the dis- 
tribution of electrostatic potential in the system. We 
then derive an explicit expression for a pressure on the 
membranes and a disjoining pressure in the gap between 
them. Our mean-field approach is verified for monovalent 
salts by molecular dynamics (MD) simulations. Simula- 
tion data fully support our theory. 

Our paper is organized as follows: In Section II some 
general consideration concerning a theoretical descrip- 
tion of an interaction between two semi-permeable mem- 
branes are presented. Here we also describe a simplified 
linearized version of the theory. Section HI contain a de- 
scription of our MD simulation approach. In Section IV 
simulation results are presented to validate the predic- 
tions of the theory. 



the study of this regime for a future work. 

We consider a hypothetical case of an infinitesimally 
thin and rigid membrane. We further assume that both 
membranes are neutral, and axis x is directed normally to 
the surfaces with a; = at the midplane of the gap. The 
membranes are located at = h/2. Our description 
thus essentially follows that of the classical non-linear 
Poisson-Boltzmann theory, except for the fact that there 
is no charge per se on the membranes and the membrane 
surface potential builds up self-consistently: accordingly 
the distribution of the charged species is a consequence 
of the semipermeable character of the membrane, which 
leads to a steric charge separation. 



A. Non-linear theory 

1. Potential 

We first introduce the dimensionless electrostatic po- 
tentials 



(1) 



with the index {i,o} standing for "in" < h/2) and 
"out" > h/2) of the confined slab. 

The non-linear Poisson-Boltzmann (NLPB) equation 
then reads 
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(2) 
(3) 



II. THEORY 

The geometry of the system under consideration is 
shown in Fig. [TJ A semipermeable membrane is in 
contact with a solution of polyelectrolyte composed of 
cations with an effective charge Z and concentration 
C, anions with charge z and concentration c. We as- 
sume here that the polyelectrolyte (here cations) cannot 
permeate through the semi-permeable membrane, while 
their counter-ions (here anions) are free to pass through 
it. The membrane is at distance h from another mem- 
brane. 

To make the formulas as transparent as possible we 
keep our analysis at the mean-field level by using the 
Poisson-Boltzmann approach. This means that we treat 
ions as point-like and neglect their correlations. In partic- 
ular, while the results of the Poisson-Boltzmann theory 
to be discussed below can be computed for any valence Z 
of the macromolecules, correlations between macro-ions 
should be taken into account in the limit of large charges 
Z in order to obtain quantitative predictions. However, 
based on earlier results [l^ , one does not expect the main 
physical picture to be altered in this case, and we leave 



where the inner inverse Debye screening length, k^, is 
defined as k\ = AttIbCoc with £b = z^e^ / [A-Kee^kBT) 
the Bjerrum length, Z ~ Z/z {< Q) is the valence ra- 
tio of large and small ions, and Coo is the concentra- 
tion of small ions far from the membrane. The outer 
inverse Debye screening length, t^Q , Cclll be defined as 
k1 = 47r^B(^^C'oo -|- Coo), where Coo is the concentra- 
tion of large ions far from the membrane. Obviously, 
it represents the inverse Debye length of the bulk elec- 
trolyte solution. Since the electroneutrality condition 

ZCao + zcao = is employed. We stress, 

however, that for this particular problem, the main ref- 
erence length scale that determines the behavior of the 
system is and this is reflected in the analysis below. 
The first integration gives a differential equation for (po 

where the integration constant is determined by the 
boundary conditions at infinity: (/)o — > and dx4>o — >■ 0. 
Imposing them on Q yields Aq = — Thus, the outer 
solution (j)o can be obtained in terms of the membrane 
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"surface" potential (ps = (j){h/2) 



2(exp[-(p\ - 1 



i(exp(-Z<^)-l)) 



2 



(5) 



For the inner compartment, (jji, the first integration of 
the NLPB equation yields 



2k? I dx 



(6) 



Here we used dx4>{x = 0) = implied by symmetry. This 
leads to 



0™ ^2 (exp[-0] - exp[-0,„]) 



(7) 



with (j)m, the (dimensionless) potential at the center of the 
film between membranes. This integral can be evaluated 
exactly as 



2 arctan 





4>m 


exp 
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leading to the Gouy-type expression 



(8) 



This value represents the bulk Donnan potential. 

Similarly, the asymptotic behavior for can be ob- 
tained from Eq. (fTTj) . Since 4>s is bounded by a con- 
stant, the condition 4>m 3> 1 as k^/i 3> 1 imposes that 
\^g-0„/2^^| ~ 7r/2. This leads to 
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(14) 



In the thin gap limit, Kih ^ 1, where the inner ionic 
clouds strongly overlap, both (j)^ and (/)s vanish. Such a 
situation would be realistic for very dilute polyelectrolyte 
solutions and/or very thin gap. One can easily verify that 

(j)m — 4>s <X Kih. 

Thus, the convergence of two semipermeable mem- 
branes is necessarily accompanied by the decrease in 
the absolute value of their potential. The idea that the 
constant potential condition is not appropriate for fully 
permeable charged membranes has been suggested be- 
fore [l^, [2l| . Now we have shown that the potential of 
neutral semipermeable surfaces should inevitably change 
and can even vanish as a result of their approach. Such a 
finding might be especially important for biomembranes, 
where an alteration of a surface potential can lead to a 
characteristic biological response. 



2. Osmotic and disjoining pressure 



(l)i{x) = -f In 



(9) 



The continuity of the electric field at the membrane 
surface, Eqs. (H])-® (no surface charge) leads to 
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Z 



(10) 



Altogether, the membrane potential (f>s and mid-plane 
potential (j)^ are thus given by the self-consistent equa- 
tions 



In 



(11) 



and 



The force balance in each part of the membrane (in 
and out) can be written 



Vp + pcE = 



(15) 



with pc the charge density and E = —dx(f> the local elec- 
tric field. Using the Boltzmann expressions for the charge 
densities in terms of the local electrostatic potentials al- 
lows to integrate this equation once. 

In midspace between the membranes (|a;| < h/2), this 
leads to 



(16) 



with c{x) = Coo Gxp[— (/)(a:)] the counter-ion concentration 
and pq a constant. 

In the outer space (|a::| < h/2) one gets 



= - In 



1 + I e 
Z 



- 1 



(12) 



In the general case, the derived equations should be 
solved numerically while the in the limits of large and 
small Kih we can also find the asymptotic analytical ex- 
pressions. In the thick gap limit, Kih ^ 1, the midplane 
potential diverges, </>,„ — ?> oo, and the equation for ipg, 
Eq. ([T^ . can be simplified to give 



Poix) = kBTc{x) + ksTC^x) + pl 



(17) 



with C{x) = Coo exp[—Z(j>{x)] the concentration of large 
ions and pl the pressure of pure solvent. 

At the membrane, there is a pressure drop due to the 
repulsion force acting on the polyelectrolyte by the mem- 
brane (and proportional to the difference of polyelec- 
trolyte concentration on the two sides of the membrane), 
i.e.: 



-iln(l-Z) 
Z 



(13) 



Pol 2 



(18) 
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This imposes po = pl, i-e. the solvent pressure, as ex- 
pected. 

The force acting on the membrane (osmotic pressure) 
can be found from the Maxwell tensor T = (P + — 
eE^E. Using V • T = 0, we find the force per unit 
surface on the membrane as Ap = T{x = 0) — T(a; = oo), 

Ap = T{x = 0) - T{x = oo) 

= keTcooil - e-'^'-) + ksTCoo (19) 

Note that by using Eq. one can demonstrate that ex- 
pression given in Eq. (|19|) is fully equivalent to Eq. (|T8|) . 
Therefore, since the disjoining pressure, 11, is defined via 



Ap = kBT{coo + Coc)-Tl, 



(20) 



it can be expressed through the concentration (or poten- 
tial) at the mid-plane as 



n = kBTCjn = kBTCooS 



(21) 



In other words, the whole effect can be expressed through 
the osmotic pressure of small ions in the mid-plane of the 
gap where the electric field vanishes. Note that the sim- 
ilar physical interpretation of the disjoining pressure be- 
tween flat solid surfaces was given long ago in the famous 
work by Langmuir Note that the disjoining pressure 
is always positive indicating an electrostatic repulsion be- 
tween neutral semipermeable membranes separating the 
similar electrolyte solutions. 

By using the expression of the disjoining pressure, (PT|) , 
we get in thick gap limit 



n ~ fcsTcoo 



2n' 



(22) 



i.e. n cx h^^. This is very similar to the famous Lang- 
muir result, but here for an a priori uncharged semi- 
permeable membrane. 
In contrast, at small Kjh 



n « knTca 



(23) 



This suggests that the osmotic pressure in the gap is 
dramatically reduced compared to the value expected in 
the bulk: 



Ap = ksTCoo = 



1 



1 



(24) 



with = ksT^Coo + Coo) the "ideal" (bulk) osmotic 
pressure. Accordingly, 11 ~ — Z/(l — Z)pi^ in this limit. 
This is one of the key results of our work. 



B. Linearized theory 

At low charge densities and low values of the electric 
potential, the description of the problem can be simpli- 
fied by linearization of the Poisson-Boltzmann approach 



(LPB). The linear approximation for the local concen- 
trations reads: Ci^o{x) = Coo (1 — ^i,o(a^)) , and Co{x) = 
Coo (^1 — Z(/)o{x)j . Substituting them into Eq. |31 we get 
its linearized version with the straightforward solutions 



! exp [ko (h/2 



1) 



cosh(Kia;) 
cosh {Kih/2) ' 



(25) 



(26) 



The electroneutrality of the membrane allows us to de- 
duce 



(27) 



Ki + Ko coth (Kih/2) 



Ko cosh (Kih/2) + Ki sinh (Kih/2) 



(28) 



In the thin gap limit, the membrane potential vanishes 
similarly to the non-linear case. For large gaps, we get 



Ki 



1 



i + vT 



(29) 



The value of this bulk Donnan potential is different from 
predicted by NLPB theory, Eq. (|T3)) . 

We also note that the potential in the midplane at large 
Kih saturates, asymptotically approaching 1 (in contrast 
to its divergence in the NLPB theory, cf. Eq.([T3])). This 
is an evidence of a failure of the linear theory in cal- 
culation of potentials (and relevant ion profiles in the 
system). 

Motivated by recent analysis [l^,[2^[2^ we then obtain 
the following expression for the pressure 



Po^PL + kBT(Coo + Coo) + 



Pi = PL+kBTCm-ZeCoo(4>i-(f>7n) + 



(30) 



(31) 



which contains quadratic terms to provide thermody- 
namic self-consistency of the linear theory. 

Using these expressions, one can calculate the pressure 
difference on both sides of the membrane to obtain [251 



Ap = pid - kBTc„ 



(32) 



As one can easily see, all the electrostatic terms cancel 
and the whole effect is expressed by the osmotic pressure 
of small ions in the mid-plane of the gap. In other words, 
physically we have arrived to the same result as in the 
NLPB case. Clearly, in the LPB case Cm = Ci(x = 0) is 
very different, so that the disjoining pressure takes the 
form 



n fcsTco 



Ko cosh (Kih/2) + Ki sinh (Kih/2) 



(33) 
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For large Kih, one can note some similarity to a repul- 
sion of solids [3, , which exponentially decays to zero 
as 



n ~ fcsTco 



■ exp 



Kih 



which obviously differs from NLPB result, Eq. 
small Kih we again get Ea. ((23| . 



(34) 
For 



III. SIMULATIONS 

The Langcvin dynamics (MD) simulations were per- 
formed on the level of the primitive model with ex- 
plicit lar ge a nd small ions using the ESPResSo simulation 
package [23 ■ We constructed a ID-periodic setup with 
two membranes fixed perpendicular to the x-axis. The 
membranes were made impermeable for cations, but "in- 
visible" for anions [28j . 

For an initial illustration of our approach we here use 
a monovalent electrolyte, and ionic species were repre- 
sented by Lennard-Joncs spheres with a central charge 
Z = 1 or z ~ ~l. We used the repulsive Lennard- Jones 
(RLJ) potential with the cut-off distance rc = 



ULjir) 




(35) 



where r is the distance between centers of two particles. 
The energy parameter e controls the strength of the in- 
teraction, and its value was fixed to e = LOksT. The 
units of length and energy in all presented data were set 
by a and e, respectively (LJ units). The bead sizes were 



set as CT- 



pp 



^ pc 



1.0. 



The interaction of ions with the membrane was set by 



Ulj{x) 




(36) 



with the cut-off distance Xc ~ 2^/^(T. 

The solvent was treated as a homogeneous medium 
with a dielectric permittivity set through the Bjerrum 
length. The electrostatic interaction between the ionic 
species was modeled by the Coulomb potential 



CouU' ij 



(37) 



where g,; = ±1. In all simulations we used = 1. 

We modeled the systems in a rectangular cell with side 
L ranging from 100 to 700 for different electrolyte concen- 
trations. The number of ions in the cell was varied from 
Np = 500 to 15000, and an equivalent number of counte- 
rions was added. The number of ions in each simulation 
was chosen to keep the ionic concentration in center of 





FIG. 2: Side view of the film confined between two semi- 
permeable membranes. Cations are shown by blue spheres, 
anions by the red spheres. One can see that anions leak out 
into the gap between the membranes. 



the membrane fixed at Coo = 0.03 and to vary the value 
of Kih in a very large range, from 0.7 to 30. 

Three-dimensional periodic boundary conditions were 
used. For Coulomb interactions we used the P3M al- 
gorithm with maximum relative accuracy of 10"'^. A 
snapshot of the system is presented in Fig. [2l 

Pressure has been evaluated via integration of the LJ 
force of cations, acting on the membrane walls 



F{x) = 4£ 



12crl 



6cr^ 



(38) 



Using this force expression, we calculated the pressure as 



h/2 



C{x)F{x)dx. 



(39) 



IV. RESULTS AND DISCUSSION 

In this section we present results of MD computer sim- 
ulations and some example calculations based on the gen- 
eral NLPB theory as well as the analytical LPB results. 

The distribution of the electrostatic potentials, (pa and 
<f>m, is shown in Fig. |3] versus Kih (symbols). Also in- 
cluded arc the exact theoretical curve, calculated with 
Eas. pT|) - p^ (solid curves). The agreement is excellent 
for all Kih, even for very large values, confirming the va- 
lidity of the mean-field approach for our system. Asymp- 
totic results are in agreement with the numerical calcula- 
tions presented in Fig. [3] At K^ft, <C 1 (strong overlap of 
an inner double layer), the membrane potential vanishes. 
Such a situation would be realistic for very dilute solu- 
tions and/or very thin gap. Another asymptotic limit of 
large films and/or concentrated solutions, Kih ^ 1 (no 
overlap of the inner ionic layers) gives Eq. ([T^ , which is 
fully supported by the simulation results. 

These pressure trends are illustrated in Fig. |31 Simu- 
lations show that at large Kih the pressure, Ap/pij sup- 
ported by the membrane, is close to the osmotic pressure 
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NLPB 




FIG. 3: Midplane and surface potentials, (j>m (top curve), (f>a 
(bottom curve), as a function of Kih. The dashed hne is the 
asymptotic behavior for 4>m in the large Kih regime according 
to Eq. (|14p : (jim. ~ 2 log [Ki/i\/2/27r] . Inset shows zoom on the 
behavior of the surface potential versus Kth. The dashed line 
here is the asymptotic behavior for 0s in the large Kih regime 
according to Ea. ([T3|l : 4>a ~ —-^ log(l — Z). In these plots we 
use Z = 1. 



of the corresponding bulk solution lS.p/pid ~ 1, and we 
deal with the standard bulk Donnan equilibrium. In this 
situation the disjoining pressure is negligibly small. At 
smaller Kih the pressure exerted on the semipermeable 
wall is much less than that in the bulk, and at very small 
Kih it approaches a constant, which is equal to the bulk 
osmotic pressure of large ions, Eq. ([M|) . This is accompa- 
nied by an increase in the value of a disjoining pressure 
in the gap. The results for the pressure and the dis- 
joining pressure obtained in the NLPB theory and simu- 
lations coincide confirming the validity of the mean-field 
approach for our system. Fig.|4]also includes the theoret- 
ical curves calculated within LPB theory. The agreement 
is quite good at very small and very large Kih, but at in- 
termediate values of Kih there is some discrepancy. The 
discrepancy is always in the direction of the pressure on 
membrane is larger than than "measured" in simulations 
and predicted by the NLPB theory. Correspondingly, 
the disjoining pressure is smaller. Still, LPB and simu- 
lation pressures are in surprisingly good harmony, espe- 
cially taking into account the simplicity of the model and 
the complexity of the system. Obviously, the effects are 
somehow included in the quadratic term of the expres- 
sion for a pressure which provided a self-consistency of 
the LPB theory. These questions, however, deserve fur- 
ther investigation and remain subject of a future work. 

Finally, we note that all the results derived above hold 
also for the case of a membrane placed at a distance 
h/2 from the neutral wall suggesting that the tight ad- 
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FIG. 4: Pressure acting on a membrane and a disjoin- 
ing pressure normalized by the ideal osmotic pressure, pid = 
kBT{coo + Coo). Solid curves show NLPB predictions, dashed 
curve shows LPB theory predictions, symbols present simula- 
tion results. 



hesion of semipermeable membrane to the neutral wall 
is impossible within our scenario. In general, for more 
complex systems, where the membrane adhesion is con- 
trolled by the competition of several effects, the physical 
mechanism we considered here should reduce the attrac- 
tive interactions by orders of magnitude similarly to what 
was predicted for other types of electrostatic interactions 
in the membrane systems [1^ H^l- It would be also 
worthwhile to emphasize that our derivation can easily 
be modified for the situation where a thin film separates 
the reservoirs with oppositely charged polyelectrolytes. 
In this case, however, no attraction between semiper- 
meable membranes takes place as it would be tempting 
to expect. The point is that only solution for such a 
configuration is = 0. Hence, all phases are homoge- 
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neous and neutral. This result can, however, still be of 
help when new synthetic delivery systems are designed. 
For example, to avoid repulsion from the cell membrane 
a semi-permeable neutral container should contain posi- 
tively charged molecules of drugs or proteins. 

To summarize, we have examined theoretically the sit- 
uation of an interaction of two neutral semipermeable 
membranes separated by a thin film. Our mechanism 
predicts an alteration of the membrane potential during 
the approach, a decrease in osmotic pressure on mem- 
branes when they are in a close proximity, and an electro- 
static repulsion between them. Our analysis also allows 
one to express a disjoining pressure in the film through 



the osmotic pressure of counter-ions in the midplane. 
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